Express Mai!: EV355345809US j 
Date of Deposit: September 19, 
2003 ! 

System and Method for Probabilistic Criticality 
Prediction of Digital Circuits 

FIELD OF THE INVENTION 

5 

This invention relates to design automation of digital integrated circuits. More 
specifically, it relates to static timing analysis of digital circuits in the presence of delay 
variations. Yet more specifically, the invention relates to determining probabilities of 
individual nodes and edges of the circuit's timing graph being critical, and to 
1 0 probabilistic enumeration of critical paths. 

RELATED APPLICATIONS 

C. Visweswariah, "System and Method for Statistical Timing Analysis of Digital 
1 5 Circuits," Docket number YOR92003-0403US 1 , U.S. Patent Application Number (to be 
assigned) filed on 09/18/03. 

C. Visweswariah, "System and Method for Incremental Statistical Timing Analysis of 
Digital Circuits," Docket number YOR92003-0401, U.S. Patent Application Number (to 
20 be assigned) filed on 09/1 8/03. 

The descriptions set forth in these co-pending applications are hereby incorporated into 
the present application by reference in their entirety. 
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BACKGROUND OF THE INVENTION 



With each succeeding generation of integrated circuit technology, variability is 
5 proportionately increasing. The sources of such variability include manufacturing 
variations, device fatigue, environmental variations and phase-locked loop (PLL) 
variations. In the case of manufacturing variations, the front-end-of-the-line (FEOL) 
which are the layers that define the active transistors show variation in the transistor's 
electrical characteristics. Physical quantities such as the length of the gate, depth of the 

10 semiconductor junction or thickness of the oxide cannot be perfectly controlled during 
manufacturing and hence show variations, which lead to variations in the behavior of the 
transistors. As the physical dimensions get smaller in modern technologies, variability is 
proportionately increasing. In addition, the back-end-of-the-line (BEOL), which consists 
of the metal interconnect layers, also exhibits variability. For example, the thickness, 

1 5 width and inter-layer dielectric thickness of each metal layer are sources of variability. 
These in turn cause the wires to change their delay, and in fact these sources of variability 
can change the delay of gates which are driving them and gates which are driven by them. 

The second main type of variations is due to device fatigue effects such as hot electron 
20 and negative bias temperature instability (NBTI). After a long period of use in the field, 
transistor characteristics change due to these physical phenomena, leading to changes in 
the delay of circuit components. 
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The third main type of variations is due to environmental effects such as temperature and 
power supply voltage. 

The fourth main type of variations is PLL variations which can include PLL jitter and 
5 duty-cycle variability. 

It is to be noted that in addition to the above, there are other sources of variation such as 
model-to-hardware miscorrelation, silicon-on-insulator (SOI) history effects and coupling 
noise. These other types of variation can also be considered during statistical timing 
10 analysis of digital integrated circuits. 

The variation of delays shown by gates and wires in an integrated circuit can be classified 
in many different ways. The variation may be from batch-to-batch during the 
manufacturing, wafer-to-wafer, chip-to-chip or within a single chip. Lens aberration 

1 5 effects during photolithography, for example, can cause variation of the effective length 
of transistors across a reticle field. There can be temperature and power supply voltage 
variations across a chip. The variations can also be classified by the time scales during 
which variability develops. For instance, fatigue effects cause variability over a period of 
years, whereas across the chip temperature or power supply gradients can develop over 

20 seconds or milliseconds, and coupling noise variations can occur in nanoseconds or 

picoseconds. Whichever way they are classified, it is abundantly clear that these sources 
of variation are making integrated circuit analysis and design more difficult and must be 
accurately accounted for during timing analysis. 
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The traditional timing methodology to handle such variability is to conduct multiple static 
timing analyses at different "cases" or "corners" to determine the spread of performance 
of the circuit under these variations. Corners may include, for example, "best case," 
5 "nominal" and "worst case." Unfortunately, the traditional methodology is breaking 
down because the number of independent and significant sources of variation is 
numerous, and too many timing runs would be required. One way to combat this is to 
worst-case or guard-band against some sources of variation, but this causes pessimism in 
the performance prediction. Another way to combat the explosion of timing runs 
10 required is to skip the analysis at certain corners, but this is risky since the performance 
of the circuit may be unacceptable at the skipped corners and this may be manifested by 
chips failing on the tester or in the field. Because of these effects, traditional timing 
methodologies are rapidly becoming burdensome, as well as risky and pessimistic at the 
same time. 

15 

A solution to the problems faced by traditional timing methodologies is statistical or 
probabilistic timing analysis. In such an analysis, timing quantities such as delays, arrival 
times and slacks are not treated as single numbers, but rather as probability distributions. 
Thus the full probability distribution of the performance of the circuit under the influence 
20 of variations is predicted by a single timing run. The problems of unnecessary risk, 

excessive timing runs and pessimism are all potentially avoided. Four examples of such 
statistical timing methods in the prior art include Liou et al [J-J. Liou, K-T. Cheng, S. 
Kundu and A. Krstic, "Fast statistical timing analysis by probabilistic event propagation," 
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Proc. Design Automation Conference, June 2001, Las Vegas, NV, pages 661 — 666], 
Scheffer [L. Scheffer, "Explicit computation of performance as a function of process 
variation," Proc. ACM/IEEE workshop on timing issues in the specification and synthesis 
of digital systems, December 2002, Monterey, CA, pages 1 — 8], Gattiker et al 

5 [A. Gattiker, S. Nassif, R. Dinakar and C. Long, "Timing yield estimation from static 
timing analysis," Proc. IEEE International Symposium on Quality Electronic Design 
(ISQED), 2001, pages 437—442] and Jess et al [J. A. G. Jess, K. Kalafala, S. R. Naidu, 
R. H. J. M. Otten and C. Visweswariah, "Statistical timing for parametric yield prediction 
of digital integrated circuits," Proc. Design Automation Conference, June 2003, Anaheim, 

10 CA, pages 932 — 937]. The references cited above are herein incorporated by reference 
in their entirety. 



PROBLEMS WITH THE PRIOR ART 



15 There are several reasons why prior-art methods are not suitable in practice. The main 
shortcoming of prior-art methods is their failure to correctly account for correlations. 
The delays of gates and wires on an integrated circuit exhibit strong correlation. 
Consider a simple example to demonstrate the importance of correlations. A chip that 
has 50,000 latches has 50,000 setup timing tests and 50,000 hold timing tests. Assume 

20 that the probability of each of these tests being successfully met is 99.99%. If the 
100,000 probabilities are perfectly correlated, then the overall probability of the chip 
working correctly is 99.99%. That is, if one test passes, they will all pass. However, if 
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the probabilities are independent, then the probability of making a working chip is 0.9999 
raised to the 100,000 th power, which is a paltry yield of 0.005%! 

There are many sources of delay correlation, and some examples are described below. 
5 Two paths may share some of the same gates; in this case, the delays of the two paths are 
correlated. If a particular chip has stronger P-type transistors due to some manufacturing 
variations, chances are that every single P-type transistor on that chip will be stronger. If 
the 3 rd level of metal is a little thicker, it is likely to be thicker across the entire chip. A 
launching path (path that gets the data signal to a latch) and a capturing path (path that 

10 gets the corresponding clock signal to the same latch) may exhibit some commonality 
and therefore have correlated delays. The commonality could include sharing some gates 
along the path, sharing metal layers, sharing a power supply voltage island, sharing gates 
of the same type, etc. Delays of gates may also be correlated because of their physical 
position on the surface of a chip. For example, two gates that are close to each other are 

1 5 unlikely to see significant relative channel length variation and their delays are therefore 
likely to be tightly correlated. 

All of these sources of variation must be accurately taken into account, or else the results 
will not be meaningful. In Liou et al, the delay of each individual gate is considered to be 
20 completely independent of any other gate, rendering the analysis unusable in practice. In 
Gattiker et al, gate delay correlations are considered. Critical paths are enumerated, but 
when the effects of these paths are combined, the delay of each path is considered to be 
independent of any other, thus ignoring several important sources of correlation. 
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Among the prior art methods, one that takes into account correlation is Jess et al. This 
method is a path-based method. The delay or slack of each path is collected and 
represented as a first-order model of the sources of variation. Then the slack of the 
5 overall circuit or chip is computed by combining these path slacks in a correlated 
probabilistic fashion. Unfortunately, there are an exponential number of paths in any 
circuit. It is not realistic to list and analyze all the paths in the circuit. Jess et al suggest 
that the top N critical paths be considered, but of course there is no guarantee that the 
(N + \) sl path (or any path other than the first N paths) will not be critical at some point 
1 0 or comer in the process space or space of variations. All path-based methods have the 
fundamental limitation that the number of paths is too large and some heuristic must be 
used to limit the number of paths submitted for detailed analysis. 

It is to be noted that although there are many significant sources of correlation in the 
1 5 delay variability of integrated circuits, there are some completely random sources of 

variation as well. For example, the oxide thickness in transistors in a modern technology 
is only a few atoms thick and for various reasons it is possible for transistors to have one 
more or one less layer of atoms, leading to variations that are quite random from 
transistor-to-transistor. While the prior-art method of Liou et al can handle such random 
20 variations, other prior-art methods such as those of Gattiker et al and Jess et al cannot. 

Further, all of the above prior-art methods have a number of other shortcomings. As 
soon as timing analysis is completed, the next natural question is, "What is the critical 
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path?" so that optimization effort to improve the circuit's timing performance can be 
appropriately directed. None of the prior-art methods describes how to obtain an ordered 
list of the critical paths, nor do they predict the probability of a path being critical. They 
do not obtain the probability of any node or edge of the timing graph of the circuit being 
5 critical. 

For these and other reasons, statistical timing methods that have been proposed in the 
literature are not used in industrial practice. 

10 ASPECTS OF THE INVENTION 

An aspect of this invention is a system and method to determine the criticality probability 
of each edge, node and path of the timing graph of a digital circuit, in the presence of 
delay variations. 

15 

Another aspect of this invention is a method for criticality probability determination with 
linear complexity in the size of the timing graph and the number of sources of variation. 

Another aspect of this invention is a method that determines the single most likely critical 
20 path of a digital circuit. 

Another aspect of this invention is a method that enumerates the most likely critical paths 
of a digital circuit along with their probabilities of being critical. 
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Another aspect of this invention is a method for criticality probability determination that 
takes into account correlations between delays of individual gates and wires, and 
correlations between delays of paths of the circuit. 

5 

Another aspect of this invention is a method for criticality probability determination that 
allows delay models that contain a deterministic part, a correlated random part and an 
independent random part. 

10 

SUMMARY OF THE INVENTION 



The present invention is a system and method for determining criticality probability of 
each node, edge and path of the timing graph of a digital circuit in the presence of delay 

1 5 variations. The delay of each gate or wire is assumed to consist of a nominal portion, a 
correlated random portion that is parameterized by each of the sources of variation and an 
independent random portion. Correlations are taken into account. Both early mode and 
late mode timing are included; both combinational and sequential circuits are handled; 
static CMOS as well as dynamic logic families are accommodated. The criticality 

20 determination complexity is linear in the size of the graph and the number of sources of 
variation. The invention includes a method for efficiently enumerating the critical path(s) 
that is/are most likely to be critical. 
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BRIEF DESCRIPTION OF THE FIGURES 



The foregoing and other objects, aspects, and advantages will be better understood from 
the following non-limiting detailed description of preferred embodiments of the invention 
with reference to the drawings that include the following: 

Figure 1 is a block diagram of one preferred embodiment of the invention depicting 
criticality determination. 

Figure 2 is a flow chart of the preferred method of determining criticality of every node, 
edge and path of the timing graph of a digital circuit. 

Figure 3 shows a sample combinational circuit to demonstrate criticality analysis. 

Figure 4 shows the timing graph corresponding to the circuit of Figure 3, with various 
forward-propagation arrival tightness probabilities annotated on the graph. 

Figure 5 shows the timing graph corresponding to the circuit of Figure 3, with various 
node criticality probabilities and edge criticality probabilities annotated on the graph. 

Figure 6 shows the timing graph corresponding to the circuit of Figure 3, with various 
backward-propagation required arrival tightness probabilities annotated on the graph. 
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Figure 7 is a block diagram of one preferred embodiment of an output report. 



DETAILED DESCRIPTION OF THE INVENTION 

5 

An inventive statistical criticality determination flow is shown in Figure 1 . The first 
input is the netlist representing the structure of the circuit to be analyzed, shown in box 
100. The second input is a set of timing assertions, box 110. These typically include 
arrival times at the primary inputs, required arrival times at the primary outputs, 

10 information about the phases of the clock, and details of external loads that are driven by 
the primary outputs. The assertions can be in the form of deterministic numbers or 
independent probability distributions or correlated probability distributions. The third 
input is a set of parameterized delay models, box 120. These allow the timer to 
determine the delay of all components of the electrical circuits (such as gates or wires) as 

1 5 a function not only of traditional delay-model variables (like input slew or rise/fall time, 
and output load) but also as a function of the sources of delay variation. For example, a 

first-order linear model may be employed, like the one shown below: 

it 

delay = a 0 + J^.Ax, + a n+l A/?, 

m 

where the delay consists of a deterministic (constant) portion a 0 , a correlated (or global) 

n 

20 portion ^^Ajc, and an independent random (or local) portion a n+] AR . The number of 

sources of variation is n , and a n i = 1,A ,n are the sensitivities of the delay to the 
sources of variation x n i = 1,A , n , and a l}+] is the sensitivity to an independent random 
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source of variation R . The notation Ax,, denotes the deviation of x i from its mean or 

nominal value, and AR denotes the deviation of R from its mean or nominal value. It is 
to be understood that the delay models can be stored in a pre-characterization step or 
calculated on the fly as required. The format in which they are stored could include 
5 analytical delay equations or table models. The next input is information about the 
statistics of the sources of variation, box 130. This input typically contains a list of the 
sources of variation with a mean value and standard deviation for each source of variation. 
Any correlations between the sources of variation are specified here. 



10 The statistical criticality determination program, box 140, accepts all of these inputs and 
produces two types of criticality reports. Box 150 shows a report that includes the 
criticality probability of each edge, node and path of the timing graph. Criticality 
information may also contain the probability of passing each timing test at each latch or 
dynamic circuit, and the probability that each end point determines the slack of the 

1 5 overall circuit. Box 1 60 shows a report that lists critical paths in order of probability of 
being critical. Either a specified number of critical paths are reported, or critical paths 
are listed till their cumulative probability exceeds a user-specified threshold. It is to be 
understood that the circuit being timed can be very large and can consist of millions of 
gates and wires. All the information mentioned above could be extremely voluminous, 

20 and so it is ordinary practice to provide options to selectively output the required 

information, or even calculate or graphically display all this information on-demand. 
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Next, the inventive criticality determination (box 140 of Figure 1) is explained in detail in 
reference to Figure 2. The first step (box 210) is to read the netlist that contains details of 
the topology of the circuit to be timed. Since this netlist is often in hierarchical form, it is 
flattened (i.e., the number of levels of hierarchy is reduced). The assertions (each of 
5 which can be either deterministic or probabilistic) are read and so are the parameterized 
delay models. Information about the sources of variation such as the mean and standard 
deviation values and any correlations between the sources of variation are inputted and 
stored in memory in box 210. 

1 0 The next major step, shown in box 220, is the construction of the timing graph, which is a 
step that is familiar in all static timing programs. In the graph, each node represents a 
node or signal of the circuit and each arc or edge represents a delay in the circuit incurred 
when a logical transition (from low to high or high to low) is transmitted through a circuit 
component such as a gate or wire. All possible valid logical transitions of the circuit are 

1 5 therefore captured in this graph. Arrival times are typically stored on the nodes of the 
graph and delays of individual gates and wires on the edges of the graph. An arrival time 
in late mode is the earliest time at which the corresponding signal is guaranteed to be 
stable at its correct logical value, and an arrival time in early mode is the earliest time at 
which the corresponding signal can change from its previous cycle stable logical value 

20 (i.e., the output cannot change later than the late mode arrival time and cannot change 

earlier than the early mode arrival time). Sequential elements and dynamic circuits in the 
graph are represented by a special kind of edge called a test segment, which is an 
indication to the timing program that a timing test must be performed between two nodes 
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of the graph to ensure correct timing operation of the circuit. Building of such a graph 
for both gate-level netlists and transistor-level netlists is known in the prior art. 

The next step is to conduct statistical timing analysis of the electrical circuit, box 230, 
5 which is preferably, but not necessarily, conducted as taught in the co-pending 

application: C. Visweswariah, "System and Method for Statistical Timing Analysis of 
Digital Circuits," Docket number YOR92003-0403US 1 , U.S. Patent Application Number 
(to be assigned) on 09/1 8/03. During the statistical timing analysis, arrival and required 
arrival tightness probabilities on each edge of the timing graph are saved for future use. 

10 

Next, in box 240, the tightness probabilities of all end points are determined by an 
inventive method that will be described in detail later. This is an optional step. If overall 
criticality probabilities are desired on a per-end-point basis, this step can be skipped. 
However, if overall criticality probabilities are desired for the entire electrical circuit, 

1 5 then this step is required. Next, in box 250, a backward propagation is carried out to 
determine the criticality probabilities of each edge and node of the timing graph by an 
inventive method that will be explained later. Finally, in box 260, timing reports are 
produced of the type in boxes 150 and 160 of Figure 1, either in a file or as return values 
of a function call. These reported values can include probabilities of individual timing 

20 tests being met; probabilities of primary outputs meeting their required arrival times; 
arrival tightness probabilities, required arrival tightness probabilities; node criticality 
probabilities; edge criticality probabilities; path criticality probabilities; a list of paths 
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ordered by path criticality. A list of paths in order of criticality probability is useful to 
the user or calling program in order to direct effort at optimizing the timing of the circuit. 

It is to be understood that the forward propagation phase of the regular statistical timing 
5 of box 230 can optionally be combined with the end point criticality determination step of 
box 240 for efficiency purposes; likewise, the backward propagation phase of the regular 
statistical timing of box 230 can optionally be combined with the backward propagation 
required to determine edge and node criticality probabilities in box 250 for efficiency 
purposes. 

10 

The details of the inventive calculations involved in boxes 240 and 250 are explained in 
the following paragraphs, in relation to Figures 3, 4, 5 and 6. The explanation below 
pertains only to late mode analysis, but is easily extended to early mode analysis as well. 
With each edge of the timing graph is associated an arrival tightness probability 
1 5 (determined during forward propagation of regular statistical timing analysis), a required 
arrival tightness probability (determined during backward propagation of regular 
statistical timing analysis) and an edge criticality probability. Associated with each node 
of the timing graph is a node criticality probability. The procedures below describe how 
the edge and node criticality probabilities are determined. 

20 

Figure 3 shows a simple combinational circuit to demonstrate statistical criticality 
determination. The circuit has four primary inputs labeled 301, 302, 303 and 304, one 
internal timing point labeled 305, and two primary outputs labeled 306 and 307. Figure 4 
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shows the directed acyclic timing graph corresponding to the circuit of Figure 3. The 
arrival tightness probabilities of edges A, B, C, D, E and F that were determined during 
box 230 of Figure 2 are annotated on Figure 4. The arrival tightness probabilities shown 
in Figure 4 were arbitrarily chosen and are for demonstration purposes only. Once the 
5 primary outputs 306 and 307 are reached, the forward propagation of arrival times is 
complete. We would like to determine path criticality probability as a product of arrival 
tightness probabilities along each path. So we now need to determine the probability of 
the end points of the timing graph (the two primary outputs 306 and 307 in our simple 
example) being critical, which is the task of box 240 of Figure 2. To do this, the timing 

1 0 graph is augmented with two dotted-line edges, G and H, from the respective end points 
to a fictitious graph node called the "sink" which has been labeled S in Figure 4. The 
delay of each such edge is the negative of the asserted required arrival time of the 
corresponding primary output. Now a statistical "max" operation is performed at the sink 
node of the graph for late mode analysis. Arrival tightness probabilities are computed in 

1 5 the same manner as during regular statistical timing analysis. For illustration, tightness 
probabilities of 0.5 and 0.5 were chosen. 

In this augmented graph, if a path is to be critical, the edges along the path must 
determine the arrival time of the sink node of every single edge along the path. Therefore, 
20 the criticality of each path is simply the product of the arrival tightness probabilities 

along the path. For example, path 303-C-305-E-306-G-S has a criticality probability of 
0.3x0.6x0.5=0.09 or 9%. The criticality probability can thus be computed for every path 
in this graph, and it can be observed that the criticality probabilities of all the paths add 
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up to 1 .0 or 100%. It is to be understood that this procedure can be readily extended to 
sequential and dynamic circuits where there is a dotted line edge of the appropriate delay 
introduced from each test segment of the graph to the sink node. The appropriate delay 
will be in terms of the clock arrival time and the guard time of the timing test. 

If this procedure (box 240 of Figure 2) is skipped, then each end point is assigned a 
criticality probability of 1 .0 and criticality computations continue, but on a per-endpoint 
basis. 

Next, we will focus on the procedure for box 250 of Figure 2. Referring to Figure 5, the 
backward propagation of edge and node criticality probabilities begins by noting that the 
sink node has a criticality probability of 1 .0 since all paths must pass through the sink 
node. Node probabilities are depicted in Figure 5 by ovals containing the node number. 
Since edges G and H have a 0.5 arrival tightness probability, the criticality probability of 
these two edges is also 0.5, since whichever is "tight" will be on the critical path. 
Moving on, the criticality probability of nodes 306 and 307 are simply the criticality 
probabilities of edges G and H, respectively, since nodes 306 and 307 have just one 
fanout each. Next, the criticality probability of edge E is the product of its arrival 
tightness probability and the criticality probability of node 6, i.e., 0.6x0.5=0.3. Similarly, 
the criticality probability of edge A is 0.4x0.5=0.2, and that of edge F is 1.0x0.5=0.5. 
Thus, in general, the criticality probability of an edge is the product of its arrival tightness 
probability and the criticality probability of its destination node. The criticality 
probability of node 305 is the sum of the criticality probabilities of edges E and F, i.e., 
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0.5+0.3=0.8. In general, the criticality probability of a node is the sum of the criticality 
probabilities of its fanout edges. All the node and edge criticality probabilities have been 
shown for the sample circuit in Figure 5. In general, node and edge criticality 
probabilities are determined by these simple methods during a levelized backward 
5 traversal of the timing graph. 

The backward propagation of required arrival tightness probabilities is shown on the 
same sample timing graph in Figure 6. At each node of the timing graph, the required 
arrival tightness probability (RATP) of a fanout edge is defined as the probability that 

1 0 that edge determines the required arrival time of the node. In Figure 6, each edge has 
three probabilities annotated on it. The one near the source of the edge is the required 
arrival tightness probability; the one near the middle of the edge is the edge criticality 
probability; and the one near the sink of the edge is the arrival tightness probability. 
When there is only one fanout edge, as in the case of node 306, the required arrival 

1 5 tightness probability of that edge is simply 1 .0, i.e., it determines the required arrival time 
of that node with 100% probability. At node 305, on the other hand, there is a 3/8 
probability that the required arrival time of node 305 is determined by edge E and a 5/8 
probability that it is determined by edge F. In general, the required arrival tightness 
probabilities of fanout edges of a node are in the proportion of their tightness 

20 probabilities, and must sum up to 1 .0. Once the required arrival tightness probabilities 
have been propagated to the primary inputs, a procedure analogous to that followed 
during forward propagation is followed. A dotted line edge is added from a source node 
(labeled "So" in Figure 6) to each primary input. The source node has a criticality 
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probability of 1 .0, and each of the dotted line edges has an arrival tightness probability of 
1 .0. Thus the edge criticality probability of each of these dotted line edges is equal to the 
criticality probability of the corresponding primary input, and the required arrival 
tightness probability of each edge is in turn equal to the dotted edge criticality probability. 

If all of the propagations of arrival times, required arrival times, tightness probabilities 
and criticality probabilities are performed as above, there are a number of useful 
properties of these quantities that can be exploited heavily during analysis and 
optimization of a circuit: 

1 . Each of the source and sink nodes has a criticality probability of 1 00%. 

2. The sum of the arrival tightness probabilities of all edges incident on any node of 
the timing graph is 100%. 

3. The sum of the required arrival tightness probabilities of all edges emanating from 
any node of the timing graph is 100%. 

4. The criticality probability of any edge is the product of its arrival tightness 
probability and the criticality probability of its sink node. 

5. The criticality probability of any edge is the product of its required arrival 
tightness probability and the criticality probability of its source node. 

6. The criticality probability of any node is the sum of the criticality probabilities of 
all edges emanating from that node. 

7. The criticality probability of any node is the sum of the criticality probabilities of 
all edges incident on that node. 

8. The sum of the criticality probabilities of all the primary inputs is 100%. 
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9. The sum of the criticality probabilities of all the end points (primary outputs for a 
purely combinational circuit) is 100%. 

10. The sum of the edge criticalities of the edges of any cutset of the timing graph that 
separates the source and sink nodes into separate partitions is 100%. 

5 11. The product of arrival tightness probabilities along any path of the timing graph 

is equal to the product of the required arrival tightness probabilities along the 
same path; both of these products represent the criticality probability of the path. 

12. The sum of path criticality probabilities over all paths of the timing graph is 100%. 

13. The arrival tightness probability of an edge is an upper bound on the criticality 
1 0 probability of any path that passes through that edge. 

14. The required arrival tightness probability of an edge is an upper bound on the 
criticality probability of any path that passes through that edge. 

15. The criticality probability of any edge is an upper bound on the criticality 
probability of any path passing through that edge. 

15 16. The criticality probability of any timing node is an upper bound on the criticality 

probability of any path passing through that timing node. 

17. The criticality probability of any node is the sum of the criticality probabilities of 
all paths in the fanout cone of that node. 

18. The criticality probability of any node is the sum of the criticality probabilities of 
20 all paths in the fanin cone of that node. 

19. The criticality probability of an edge is the sum of the criticality probabilities of 
all paths that pass through that edge. 
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The above properties can be exploited during analysis and optimization of a circuit. For 
example, during path enumeration, various upper bounds on path criticality listed above 
can be used to prune the list of paths as the list is generated, whether in a depth-first or 
breadth-first fashion. The list of paths may then be either reported to the user or 
5 submitted for more detailed analysis. During optimization, the edges and nodes with the 
highest criticality probabilities can be targeted for improvement. Therefore it is to be 
understood that the above quantities can be exploited in many different ways in a variety 
of contexts. 

1 0 After the criticality determination procedure is complete, the results are fed back to the 
user by means of a programming interface or timing report stored in a hard disk file or 
files. Timing reports are illustrated in box 710 of Figure 7. A timing report typically 
consists of circuit information (box 720) and criticality probability information (box 730). 
The results of the inventive criticality computations can be reported in conjunction with 

1 5 statistical timing results for easy digestion by the user. Circuit information can include a 
list of gates, components and wires; or a list of paths; or a list of nodes; or a list of 
sequential elements; or a list of end points (primary outputs and timing tests); or a list of 
clock phases; or a list of edges of the timing graph. These items can be sorted and 
filtered in various ways to make the report intuitive and productive to the reader of the 

20 report. The corresponding statistical timing criticality information in the case of a node, 
edge or path, could include the node's, edge's or path's criticality probability, 
respectively. The criticality probability could be reported on a global basis or on a per- 
endpoint basis. In addition, arrival tightness probabilities and required arrival tightness 
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probabilities are available to be reported for each edge of the timing graph. Similarly, for 
a timing test, the criticality probability of the test may be reported. It is to be understood 
that each of the criticality probabilities in the above description can be one of an early- 
mode or late-mode criticality probability; one of a rising or falling criticality probability; 
5 and a criticality probability that is specific to a particular phase of a particular clock. It is 
to be further understood that once the criticality analysis is completed, these results can 
be reported in a variety of useful ways. 

It is to be understood that the detailed description of this invention was explained in the 
1 0 context of a simple combinational circuit with a simple timing graph. One of ordinary 
skill in the art will be able to extend these concepts to accommodate separate rising and 
falling arrival times and delays; sequential circuits; circuits with transparent latches; 
extensions to handle slew (rise/fall time) propagation and effects; circuits with 
probabilistic guard times on the latch setup and hold tests; circuits with probabilistic 
1 5 guard times on dynamic circuit timing tests; circuits with multiple clock phases; 

probabilistic delay models that are purely random; probabilistic delay models that are 
correlated; sources of variation that are random; sources of variation that are correlated; 
and cases in which clock cycle periods and timing assertions are probabilistic. 

20 Given this disclosure, it is apparent to one skilled in the art that the inputs received by the 
"Statistical criticality determination" process (box 140 of Figure 1) can be any input 
generally known to computer systems, including but not limited to: keyboard or mouse 
entries, disk, tape, CD-ROM, network connection, fiber optic connection, radio frequency 
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link, infra red link, etc. Further the outputs including the "Criticality of each edge, node 
and path of the timing graph" (box 150 of Figure 1) and the "List of critical path(s) with 
criticality probabilities" (box 160 of Figure 1) can take the form of any known computer 
output. These outputs include but are not limited to: printed output from a printer, 
5 images on a graphical user interface (GUI) or CRT, content on storage media (e.g., 
memory, CD-ROM, disk, diskette), files, information transmitted over a network (fiber 
optic, telephone, cable, radio frequency, infrared, etc.). 
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